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I ABSTRACT 

This  report  is  a user  manual  for  the  program  QUADRAW 
which  draws  visible-line-projection  pictures  of  objects 
bounded  by  sections  of  quadric  surfaces.  The  report 
describes  three  versions  of  the  program,  one  each  for  the 
CDC  6600,  Univac  1110  and  IBM  System/360-370  computers. 

It  also  provides  the  necessary  documentation  to  permit 
one  to  modify  the  program  for  any  other  computer. 
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1.  INTRODUCTION 
1. 1 Outline  of  the  Program 

QUADRAW  is  a computer  program  which  draws  visible-line 
pictures  of  solid  objects  bounded  by  sections  of  quadric 
surfaces.  The  algorithm  and  the  original  program  were 
developed  by  P.  Woon  [1,2]. 

QUADRAW  reads  a numerical  description  of  a scene  and 
one  or  more  vantage-point  specifications  from  an  input  file 
[Fig.  1.1].  The  program  produces  a drawing  of  the  scene 
from  each  vantage  point  [Fig.  1.2].  Drawings  are  output  on 
a digital  plotter.  Each  drawing  of  a scene  may  show  (1)  all 
lines,  (2)  hidden  lines  dashed,  and  (3)  only  the  visible 
lines.  Perspective  or  orthographic  projection  of  the  scene 
may  be  selected  at  the  option  of  the  user. 

QUADRAW  is  written  in  Fortran  and  consists  of  a main 
program  and  more  than  70  subroutines  and  functions.  At  the 
present  time,  there  are  available  versions  for  the  CDC  6600, 
Unlvac  1110  and  IBM  360/370  computers.  It  should  be  quite 
simple  to  adapt  one  of  these  three  versions  to  any  other 
computer  having  a Fortran  compiler  and  a large  memory  space.* 

The  QUADRAW  source  text  is  maintained  as  a library  file 
on  a CDC  6600.  This  file  contains  all  versions  of  the  pro- 
gram. The  CDC  text  editor  UPDATE  can  generate  the  text  of 


j 

I 


Including  system  and  plotting  support,  the  program  re- 
quires 150k  words  on  a CDC  6600,  I850O/38OOO  (I/D)  words 
on  a Unlvac  1110  and  250k  bytes  on  an  IBM  36O  computer. 


1 

2 


3 


•t 


5 

6 


S 


20.0 

STACK  QF  8 

CUT  SPHERES 

l.OE-03 

1.0E-0<. 

l.OE-05 

l.OE-09 

l.OE-02 

1 

1 

1 .000 

l.OlO 

1.000 

AU  ♦9»  ♦A.  *9  ) 

-8.000 

-0.000 

-8.000 

32.000 

2 

1 

l.COO 

l.OCO 

1.000 

ATI»9.*9»-9I 

-8.000 

-8.000 

8.000 

32.000 

3 

1 

1.000 

1.000 

1.000 

AT(99,-9.+9l 

-3.000 

8.000 

-8.000 

32.000 

1 

1.000 

1.000 

1.000 

AH+9»-9»-9» 

-8.000 

8.000 

8.000 

32.000 

5 

1 

1.000 

l.COO 

1.000 

AT(-9.+9,»9l 

8.000 

-8.000 

-0.000 

32.000 

6 

1 

1.000 

l.OOO 

1.000 

AT(-9. *9,-91 

0.000 

-8.000 

8.000 

32.000 

7 

1 

1.000 

1.000 

1.000 

AT(-9,-9,*9l 

8.000 

8.000 

-0.000 

32.000 

8 

1 

1.000 

1.000 

l.OOC 

AT(-9,-9,-9) 

3.000 

8.000 

8.000 

32.000 

9 

1 

PLANE 

1.000 

-6.000 

AT  2»6 

10 

1 

PLANE 

1.000 

2.000 

AT  Z.-2 

1 

-9 

9999 

2 

-10 

9999 

3 

-9 

9999 

' 

9 

-10 

9999 

' 

5 

-9 

9999 

6 

-10 

9999 

7 

-9 

9999 

8 

-10 

9999 

9 

-1 

9990  -3 

9990  -5 

9998  -7 

9999 

10 

-2 

9998  -9 

9998  -6 

9998  -8 

9999 

1 

9 

2 

1.000 

-9. COO 

2 

10 

2 

1.000 

-9.0u0 

3 

9 

2 

1.000 

9.000 

9 

10 

2 

l.OCO 

9.bu0 

5 

9 

e 

1.000 

-9.000 

6 

10 

2 

1.000 

-9.C0C 

7 

9 

2 

1.000 

9.000 

8 

10 

2 

1 .000 

9.000 

3 

2C.0C0 

0.000 

90.0C0 

b.OOO 

3 

2C.OOO 

0.000 

70.0C0 

30.000 

3 

2C.OOO 

-55.000 

25.000 

75.000 

3 

20.000 

10.000 

-75.000 

-55.000 

3 

20.000 

10.000 

10.000 

O.OUO 

3 

20.000 

-60.500 

-10.000 

5.  OUO 

3 

20.000 

90.000 

-50.000 

22.500 

3 

20.000 

98.000 

-22.500 

-90.000 

1 


Pig.  1.1  An  example  of  a QUADRAW  input  file  (Scene  made 
of  8 spheres  which  are  cut  by  2 planes) 
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any  implemented  version.  This  arrangement  simplifies  the 
maintenance  of  the  program  since  any  modifications  of  the 
program  will  automatically  be  applied  to  all  implemented 
versions.  Note  that  the  three  versions  differ  only  in  small 
parts  which  are  very  unlikely  to  be  changed.  A version  for 
any  other  computer  can  be  easily  added  to  this  master  source 
text . 

Although  QUADRAW  can  also  handle  objects  bounded  en- 
tirely by  planar  surfaces,  a different  algorithm  [3]  and 
program  [4],  specifically  developed  for  planar-surfaced 
objects,  would  be  faster  and  require  less  memory  space,  and 
is  recommended  in  that  case. 

1 . 2 Definition  of  Terms 

A scene  is  set  of  one  or  more  bodies . Each  body  is  a 
closed,  connected  space  in  3D,  delimited  by  a number  of 

i 

quadric  surfaces.  The  terms  object  and  body  are  used  inter-  i 

changeably . 
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A surface  bound  Is  a quadric  surface  used  to  delimit 
the  sections  of  a quadric  surface  that  bound  a body. 

A surface  intersection  is  the  locus  of  all  points  of 
two  quadric  surfaces  satisfying  Q^(x,y,z)  = Q2(x,y,z)  = 0. 

A real  vertex  is  a point  in  3-space  where  three  or 
more  quadric  surfaces  intersect.  If  these  three  or  more 
quadric  surfaces  intersect  in  several  points,  then  these 
points  are  called  a set  of  real  vertices. 

A virtual  edge  or  limb  is  the  locus  of  all  points  where 
a quadric  surface  is  tangent  to  the  lines  of  projection. 

A virtual  vertex  is  the  Intersection  in  3-space  of  a 
virtual  edge  and  a surface  intersection. 

An  edge  is  either  a virtual  edge  or  a surface  inter- 
section . 

A polar  plane  is  a plane  in  which  all  virtual  edges  of 
a quadric  surface  lie. 

A cutting  plane  is  a plane  which  intersects  the  inter- 
section of  two  quadric  surfaces. 

An  object  coordinate  system  is  a 3D  cartesian  coordin- 
ate system,  x"  y"  z" , in  which  all  object  equations  and 
vantage  points  are  specified. 


I 
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2.  DESCRIPTIONS  OF  QUADRAW  ROUTINES 


In  this  section  a brief  description  of  each  QUADRAW 
routine  is  given.  Pig.  2.1  shows  the  interconnections 
among  the  routines. 

2 . 1 Main  Program 

2.1.1  QDRAW  - controls  the  execution  of  the  pro- 
gram. It  calls  the  main  function- 
al subroutines,  opens  and  closes 
plotting  output  and  measures  exe- 
cution time  of  the  main  steps  of 
the  program. 

2.2  Main  Functional  Subroutines 


2.2.1 

INPOB 

reads  scene  specifications  and 
creates  Tables  1-5  for  the  scene. 

2.2.2 

RVERT 

computes  the  locations  of  the  real 
vertices  of  the  objects  and  stores 
them  in  Table  4. 

2.2.3 

SINT 

computes  all  surface  intersections, 
point  by  point  and  stores  them  in 
Table  6. 

2.2.4 

INPVW 

reads  vantage  point  specifications. 

2.2.5 

TRNSF 

transforms  the  surface  equations 
from  the  object  specification  co- 
ordinate system  to  the  coordinate 
system  of  the  current  vantage  point. 

2.2.6 

REDGE 

transforms  surface  intersections 
from  the  object  specification  co- 
ordinate system  to  the  coordinate 
system  of  the  current  vantage  point. 
It  also  computes  all  virtual  ver- 
tices . 

2.2.7 

VEDGE 

computes  virtual  edges  - limb  seg- 
ment s . 

2.2.8 

INSIDE 

sets  up  the  surface  minimum-maximum 
table  (Table  9)  which  is  later  used 
for  envelope  tests  in  determining 
edge  visibility.  It  also  inserts 
edge  minimum  and  maximum  coordinates 
in  Table  8. 

Fir.  2.1  Interconnections  of  QUADRAW  routines 
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2.2.9 

FINAL 

determines  the  visibility  of 

all  points  of  surface  intersections 

and  virtual  edges. 

2.2.10 

DRAW 

calls  plotting  routine  to  draw  the 
projection  of  the  scene;  the  scale 
factor  of  the  drawing  is  determined 
from  Table  9 and  the  point  coordin- 
ates are  obtained  from  Table  8. 

This  subroutine  may  draw  (1)  visible 
intersections  and  virtual  edges, 

(2)  hidden  intersections  and  virtual 
edges  as  dashed  lines,  or  (3)  all 
intersections  and  virtual  edges. 

Subroutines  2.2.1 
Subroutines  2.2.4 

-3  are  called  once  per  scene. 

-10  are  called  once  per  vantage  point. 

Other 

Funct ional 

Subroutines 

2.3.1 

CUT 

called  by  SINT  to  find  a 'cutting 
point'  if  there  are  no  real  vertices 
on  a real  intersection.  The  equa- 
tion of  a 'cutting  plane'  which 
intersects  the  intersection  is 
stored  in  Table  3. 

2.3.2 

IBMD 

an  integer  function  called  to  deter- 
mine whether  a point  is  within  a 
surface's  bounds. 

2.3.3 

IEND2 

performs  the  same  operations  as  I3ND 
on  surfaces  in  the  vantage-point 
coordinate  system. 

2.3.4 

ROTPT 

called  by  REDSE  to  rotate  (and  trans- 
late, for  perspective  projections) 
each  surface  intersection  point  to 
get  its  coordinates  in  the  vantage- 
point  coordinate  system. 

2.3.5 

DIRECT 

called  by  VEDGE  to  adjust  the  direc- 
tion of  propagation  of  an  edge. 

2.3.6 

DELEV 

called  by  SINT  and  VEDGE  to  delete 
a pair  of  vertices  from  a vertex 
table . 

2.3.7 

FSTPT 

called  by  VEDGE  to  determine  the 
first  point  on  a limb  that  has  no 
vertices . 
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2.3.8 

EXTRM 

called  by  FSTPT  to  find  the  extrema 
of  a limb. 

2.3.9 

INVIS 

called  by  FINAL  to  check  whether  a 
face  of  a surface  hides  a particular 
point  from  the  vantage  point. 

2.3.10 

WDUMP 

lists  the  current  contents  of  Tables 
1-9  and  15.  It  is  called  from 
many  subroutines  whenever  an  error 
is  detected  by  the  program  in  either 
the  description  of  an  object  or  com- 
putations. 

2.3.11 

CLOSER 

lists  all  input  cards  that  were  not 
read  after  an  error  in  input  was 
found  by  INPOB. 

2.3.12 

INITO 

initializes  all  global  variables  to 
zero  before  a new  scene  is  read  by 
INPOB. 

2.3.13 

SORTl 

sorts  an  array  of  data  into  a des- 
cending sequence. 

2 . 4 Table  Management  Routines 


These 

routines 

access  the  tables  of  the  QUADRAW  data 

structure . 

They  are 

called  by  the  routines  described  in  2.2 

and  2.3.  The  small 

letters  stand  for  single  digit  numbers. 

2.4.1 

GTmn 

is  the  n’th  routine  for  getting  the 
contents  of  table  number  m.  There 
are  12  GTmn  subroutines. 

2.4.2 

GTI5I 

gets  the  contents  of  Table  15. 

2.4.3 

STOmn 

is  the  n’th  routine  for  storing  the 
contents  of  table  number  m.  There 
are  14  STOmn  subroutines. 

2.4.4 

STI5I 

stores  the  contents  of  Table  15. 

2.4.5 

STI52 

removes  the  last  entry  in  the  list 
in  Table  15. 

2.4.6 

FULL 

is  called  when  the  data  structure 
runs  out  of  space,  that  is,  when 

Table  S or  Table  9 hits  Table  15  or 
when  Table  1-9  hits  the  end  of  the 
data  structure. 

V 

} 
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2.4.7  TABLE 

2.4.8  IHTBL9 

2 . 5 Mathematical 

2.5.1  EQlll, 
EQ211 
EQ221 

2.5.2  EQ22 

2.5.3  ROOT2 

2.5.4  SUBSl 

2.5.5  NEWT 

2.5.6  PWRT  & 
POLRT 

2.5.7  ELX21 

2.5.8  ELY2 

2.5.9  DERIV 

2.5.10  DET3 

2.5.11  POLAR 

2.5.12  IPNT 


stores  elements  of  Table  9. 

is  called  by  INSIDE  to  initialize 
Table  9. 


Subroutines 


compute  the  Intersection  point (s) 

& of  three  equations,  with  '2'  stand- 

ing for  second-order  equation  and 
'1'  standing  for  first-order 
(linear)  equation. 

computes  the  intersection  points 
of  two  second-order  equations. 

solves  a quadratic  equation  of  one 
variable,  indicating  the  number  of 
roots . 

substitutes  a variable  in  an  equa- 
tion to  make  it  into  a simpler 
equation . 

computes  the  intersection  of  two 
curves  (i.e.,  the  intersection  of 
two  surfaces  in  a plane)  by  using 
the  Newton's  method. 

help  EQ22  in  solving  two  second- 
order  equations. 

eliminates  variable  x between  a 
second-order  equation  and  a first- 
order  equation. 

eliminates  variable  y between  two 
second-order  equations  of  y and  z. 
Sylvester's  method  is  used. 

computes  3^y/9z^  or  9^z/3y^  of  a 
conic  section  at  point  y,z. 

evaluates  the  determinant  of  a 
3x3  matrix. 

computes  the  coefficients  of  the 
polar  plane  of  a surface. 

computes  points  on  intersections  of 
two  surfaces.  It  is  used  to  trace 
points  on  surface  Intersections  as 
well  as  on  limbs. 
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2 . 6 Plotting  Subroutines 

QUADRAW  assumes  that  a plotting  package  Is  available  to 
draw  pictures  on  a digital  plotter  or  a similar  device.  The 
plotting  package  should  be  able  to  perform  the  following 
tasks : 

a)  open  and  close  the  plotting  output 

b)  position  the  plotter  to  a new  picture 

c)  draw  solid  lines 

d)  move  the  pen  while  It  Is  In  the  "up"  position 
The  plotting  packages  for  the  Implemented  versions  of 
QUADRAW  are  described  in  their  respective  appendices. 


2 . 7 Timing  Routine 

2.7.1  MSTIME(MS)-  is  an  integer  function  which  returns 

as  its  value  the  elapsed  program 
execution  time  in  milliseconds  since 
time  MS  (also  in  milliseconds).  This 
function  calls  a system  routine  to 
measure  the  time.  The  system 
routines  for  the  implemented  versions 
of  QUADRAW  are  described  in  their 
respective  appendices. 
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3 . DATA  STRUCTURE 

All  QUADRAW  data  structure  and  global  variables  are 
stored  in  the  blank  COMMON  block.  The  ten  tables  used  by 
the  program  are  stored  in  array  IPOOL( 32000)  which  is  also 
EQUIVALENCEd  to  POOL (32000)  and  therefore  contains  integer 
as  well  as  real  variables.  The  tables  in  IPOOL  are 
numbered  1-9  and  15-  Tables  10-14  do  not  exist  but  could 
be  added  if  the  program  were  extended, 

3 • 1 Global  Variables 

The  following  variables  and  arrays  are  contained  in  the 
blank  COMMON  block: 

PHI,  THETA,  PSI  are  the  azimuth,  elevation  and  twist  angles 
of  the  current  vantage  point. 

PSCL  is  the  picture  scale  of  the  scene  from  the  current 
vantage  point . 

IPTYP  is  the  type  of  drawing  to  be  made  from  the  current 
vantage  point  view. 

R1 , R2 , R3,  R4 , R5 , R6,  R7 , R8 , R9  are  the  nine  elements  in 

the  rotation  matrix  which  transforms  an  object  from  the  ti 

object  coordinate  system  to  the  current  vantage  point 
coordinate  system.  ' 

EPSl,  EPS2,  EPS3,  EPS4,  EPS5  are  the  comiputatlonal  tolerances  ; 

used  in  the  program. 

IHEAD  contains  an  integer  digit  from  column  one  of  the  last 
data  card  read  by  the  program.  Digits  in  the  first 
column  of  data  cards  are  used  to  delimit  various 

sections  of  data  input.  ; 

i 

1 

5 

1 

4 


-13- 


NSURFS  contains  the  number  of  surfaces  in  the  current  scene. 

TAPE  is  a logical  variable  set  to  .TRUE,  when  the  plotting 
output  is  open. 

IPBl,  IPB2,  , IPB14,  IPB15  are  pointers  to  the  15 

possible  tables  in  IPOOL.  Note  that  Tables  10  - 14 
do  net  exist . 

MAXDM  contains  the  dimension  of  the  array  IPOOL  which  is 
32000  at  the  present  time. 

YMN,  YMX,  ZMN,  ZMX  are  used  by  subroutine  DRAW  to  store  the 
maximum  and  minimum  coordinates  in  y and  z directions 
of  the  scene  in  the  projection  plane. 

SPARE(7)  is  an  array  of  seven  spare  variables. 

LISTER  is  a logical  variable  set  to  .TRUE,  if  the  input  data 
specifies  that  tables  in  IPOOL  are  to  be  listed  after 
the  current  scene  is  entered  or  after  the  current 
vantage  point  drawing  is  made. 

ABORT  is  a logical  variable  set  to  .TRUE,  when  an  error  is 
detected  during  computations  and  the  current  object 
or  vantage  point  computations  fail.  The  program  tries 
to  continue  with  the  next  scene  or  the  next  vantage 
point . 

DUMPER  is  a logical  variable  set  to  .TRUE,  when  a listing  of 
tables  in  IPOOL  is  to  be  made  by  subroutine  WDUMP  after 
the  detection  of  an  error. 

DE  is  used  for  perspective  projections  only.  DE  is  the  dis- 


tance from  the  vantage  point  to  the  origin  of  the 
object  coordinate  system. 


-1^- 


DEE  is  used  for  perspective  projections  only,  DEE  is  the 
distance  from  vantage  point  to  the  picture  projection 
plane . 

IPTNA  is  a pointer  to  the  last  location  in  IPOOL  occupied 
by  Tables  1-9. 

033CL  is  the  scene  scale. 

IPOOL(32000)  is  the  array  containing  the  tables  of  scene 

description.  This  array  is  EQUIVALENCEd  to  POOL( 32000). 

The  following  is  a description  of  all  ten  tables  in 
IPOOL-POOL.  Tables  1 to  6 are  set  up  once  per  scene,  while 
tables  7 to  9 a.nd  15  are  set  up  once  per  vantage  point. 


* The  array  contains  integer  numbers  as  well  as  real  numbers 
(floating  point).  We  mark  integer  entries  with  I’s  and 
floating  point  entries  with  R’s. 
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3 . 2 Table  1:  Surface  Equations 

There  is  one  l4-word  element  in  this  table  for  each 
surface  of  an  object. 


1 

ISTYPE 

2 

IDEG 

3 

LINK 

FACTOR 

5 

~TrD  : 

6 

A(2) 

7 

~aT31  i 

8 

1 

9 

A(5)  1 

10 

Tltl  i 

11 

A(7)  1 

12 

“OT)  i 

13 

~^J9) 

14 

A(10) 

surface  type:  l=component,  2=auxiliary 
degree  of  equation;  l=llnear,  2=quadric 
pointer  to  the  surface's  bounds  in  Table  2* 
scale  factor  = 1.0 


the  ten  coefficients  of  the  surface 
equation : 

2 2 2 

a^x  + a2y  + a^z  + a^xy  + a_xz  + 

V *7==  * V * ^9-  * ho  ■ “ 


This 

table 

is 

This 

table 

is 

written  by: 
read  by: 


STOll  and  ST012,  both  called 
by  INPOB. 

GTOll  called  by  RVERT,Si:iT 
and  TABLE, 

aT012  called  by  RVERT  and 
SINT, 

GT013  called  by  TRNSF  and 
REDGE. 


j * The  pointer  is  set  to  zero  for  self-bounding  surfaces. 

! 

1 

j 

k » - - ---- - 
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3-3  Tiable  2:  Surface  Bounds 

There  is  one  element  of  variable  length  in  this  table  for  each 
bounded  surface*  in  an  object.  Pointer  LIM  in  table  1 points  to 
this  element.  Self-bounded  surfaces  - such  as  spheres  - do  not  have 
eleinents  in  this  table.  Pointer  LH-JK  for  a self-bounded  surface  is 
set  to  zero. 


surface  numbers  of 
bounding  surfaces 

j. 

end-of-e lament 


This  table  is  written  by  INPOB.  It  is  read  by  IBND,  IBrJD2,  TRNSF  and  HEDGE. 


* Surface  bounds  are  explained  in  section  4.2. 


I 
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3.^  T^ble  3:  Surface  Intersections 

Ihere  is  one  9-word  element  for  every  surface  intersection. 


1 

2 

3 

k 

5 

6 

7 

8 

9 


lOTYP 


ISl 


IS2 


KV 


KIP 


CUTP(7) 


CUTP(8) 


CUTP(9) 


CUTP(IO) 


I 

I 

I 

I 

I 

R 

R 

R 

R 


intersection  type 

first  surface  number 

second  surface  number 

pointer  to  real  vertex  pointers  - 

pointer  to  intersection  points  - 


Table  5 
Table  6 


equation  of  cutting  plane*: 
cutp^  + cutp0y  + cutp^z  + cutp^Q=0 


Written  by; 


Read  by: 


ST031 

from 

INPOB; 

words  1-3,  6-9 

ST032 

from 

SINT  : 

word  5 

ST051 

from 

E^OB: 

word  4 

siirr, 

HEDGE 

* Indices  of  coefficients  are  set  up  for  coipatlblllty  with  quadric 
equatior.s.  Ln  section  we  refer  to  the  cutting  plane  equation 
' >-2^  ' >-0“  • 


as  c,x  + + c-,2  + c,,  = 0. 


IB- 


S'S liable  4;  Real  Vertices 

There  is  an  element  of  up  to  19  words  In  this  table  for  each  real 
vertex  in  an  object.  For  a set  of  real  vertices,  that  is  vertex  multi- 
plicity MULTIP  > 1,  the  element  is  repeated  MULTIP  number  of  times. 


vertex  multiplicity 
X,  y and  z coordinates 
of  a real  vertex 

numbers  of  surfaces  which 
form  this  vertex,  up  to 
15  surfaces  allowed 

end-of-element  marker 


Vertex  multiplicity  gives  the  number  of  real  vertices  which  are 
fomed  by  intersections  of  the  same  three  or  more  surfaces.  If 
MULTIP  > 1,  then  subsequent  elements  in  the  set  of  real  vertices  have 
MULTIP  set  to  0. 

This  table  is  written  by  ST04l  called  from  INPOB  (words  2-4)  and 
STO^P  called  fron  PVEF5T  (words  1 , 5-last ) . 


F 
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3.6  T^ble  5:  Real  Vertex  Pointers 

There  is  a variable  length  element  in  this  table  for  each  surface 
intersection  which  has  one  or  more  real  vertices  on  it.  Pointer  KV  in 
surface  intersection  Table  2 points  to  this  element.  If  there  is  not 
any  real  vertex  on  a surface  intersection  pointer  KV  is  set  to  zero. 


I 

I 

■ 

I 

999a 

I e 

pointers  to  real  vertices, 
in  table  which  lie  on 
the  surface  intersection 


9999 1 I end-of- table  5 marker 


This  table  is  witten  by  subroutine  STO51  called  by  INPOB. 


It  is  read  by  subroutine  SUIT. 


L 
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3.7  'te.ble  6:  Surface  Intersection  Points 

There  are  three  types  of  3-word  elements  in  this  table:  points, 
end-of-line  markers  and  end-of-intersection  markers.  The  first  point 
of  intersection  of  two  surfaces  is  pointed  to  by  pointer  KIP  in  Table 
3-  This  table  usually  occupies  a large  amount  of  space. 

Point: 

X,  y and  z coordinates  of  an  inter- 
section point. 

End-of-llne  marker: 

marker  indicating  the  end  of  a line 
within  a surface  intersection. 

I End-of-intersection  marker: 

marker  indicating  the  end  of  the 
last  line  of  a surface  intersection. 


This  table  is  written  by  subroutine  STO61  from  SINT. 


3.8  Table  7:  Transformed  Surface 


There  is  a 19-word  element  in  this  table  for  each  surface  of  an 


ISTYP* 

IDEG* 

LIMC2* 

FACTOR* 

LINK15 

3(1) 

3(2) 

3(3) 

3(4) 

3(5) 


type  of  surface 
degree  of  surface 
pointer  to  Table  2 (bounds) 
scale  factor  = 1.0 

pointer  to  Table  15  (virtual  vertices) 


POLRA(7) 

POLRA(d) 

POLRA(9) 

POLRAUO) 


coefficients  of  the  transformed 
surface  equation 

coefficients  of  the  polar  plane 
equation  ** 


a)  Words  1-15  are  '/rritten  by  TRMSF  except  word  5 - LIi'IK15. 
LINK15  Is  set  Initially  to  0 by  TRNSF,  set  to  -1  by 
POLAR  If  there  is  not  any  polar  plane.  It  Is  set  to 
point  to  Table  15  by  REDIE  if  there  are  virtual  vertices. 


b)  Words  16-19  (polar  plane  coefficients)  are  written  by  POLAR. 


* Copied  from  Table  1. 

**  Indices  of  coefficients  are  set  up  for  conpatlbillty  with  quadric 
equations . 
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3.9  Table  8 : Edges 

There  are  two  types  of  elements  into  this  table.  Header  elements 
which  are  11  words  long  and  point  elements  which  are  3 words  long. 

Header  Element: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 


ITYP** 

iCLAS** 

lOHITTA** 

ISNB** 

lORMTB** 

LiriKb* 

YMAX*»* 

ZMAX*** 


I 

I 

I 

T- 

X 

I 

T 

I 

R 

R 

R 

R 


Point  Element: 


1 

2 

3 


r X* 

-yF 


R 

R 

R 


1 

2 

3 


999  998.0* 
999  998.0* 
999  998.0* 


R 

p 

R 


surface  intersection  type,  0 for  virtual  edges 


edge  class:  H,,Hp,H^  or  H^ 
number  of  first  surface 
orientation:  +1= front,  -l=back 
number  of  second  surface 
orientation:  +l=front,  -l=back 
pointer  to  the  next  edge  header 


maximum  and  minimum  y and  z 
envelope 


i 


X,  y and  z coordinates  of  a point 
If  the  point  is  invisible  X is  changed  to: 
999  997.0  for  orthogonal  projection 

-X  for  perspective  projection 

by  FIMAL. 


end-of-line  marker 


1 

2 

3 


999  999.0* 
999  999.0* 
999  999.0* 


R 

R 

R 


end-of-surface  marker 


» Set  by  REDGE  and  ITEGE. 

**  Set  by  ST081  from  REDGE  and  VEDGE. 
***  Set  by  ST085  ftxxn  DBIDE. 


'1 
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Surface-Enclosing  Envelopes 
contains  two  types  of  elements.  There  is  one 
for  each  surface  of  the  object.  A primary 
words  long.  Secondary  ele.ments  are  added  as 
A secondary  element  is  10  words  long.  The  table 


3.10  Table  9: 

This  table 
primary  element 
ele.ment  is 
required . 
contains  the  coordinates  of  the 
and  pointers  to  its  edges. 

Primary  element ; 


enclosing  envelope  of  a surface 


coordinates  of  the  enclosing 
envelope 


pointers  to  edges; 
last  pointer  followed  by 
a zero  pointer 


Secondary  element: 


pointers  to  edges; 
last  pointer  followed  by 
a zero  pointer 


This  table  is  initialized  by  INTBL9  called  by  INSIDE. 

It  is  written  by  TABLE  called  from  INSIDE,  and  read 
by  FINAL  and  DR.AW. 


* A negative  pointer  points  to  the  next  secondary  element. 
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3-11  Table  15:  Virtual  Vertices 

There  is  one  4-word  element  in  this  table  for  each 


The  x,y  and  z coordinates  of 
a virtual  vertex  written  by 
HEDGE;  read  by  VEDGE. 

Pointer  to  the  next  virtual  vertex 
element  in  a list.  Last  pointer  in 
list  is  set  to  zero. 


virtual  vertex. 


1 

2 

3 

4 


iPT'rv 


D 

P 

R 

I 


Note 

This  table  is  written  backwords  at  the  same  time  that 
Table  8 is  written  forwards.  If  Table  8 or  Table  9 should 
hit  Table  15  subroutine  FULL  is  called  to  stop  further 
execution  of  the  program. 


J 
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'4.  SCENE  AND  VANTAGE  POINT  SPECIFICATIONS 

In  this  section  we  describe  the  format  of  data  that 
specifies  3UADRAV/  scenes  and  vantage  points.  Each  scene  entered 
into  QUADRAW  is  specified  by  these  five  parts: 

1.  Scene  scale,  title  and  computational  tolerances. 

2.  Quadric  equation  of  each  surface. 

3.  Bounds  of  each  bounded  surface. 

All  surface  intersections. 

5.  All  real  vertices. 

Scene  specifications  are  usually  followed  by  one  or  more 
vantage  point  specifications.  Each  vantage  point  is  specified 
by  its  position,  type  of  projection  and  type  of  drawing  to  be 
made  by  the  program.  QUADRAW  can  use  either  orthographic  or 
perspective  projection. 

The  following  pages  describe  the  format  of  scene  and 
vantage  point  data  in  detail.  80-column  punch  cards  are  used 
to  Illustrate  data  formats.  They  represent  any  30-or-more- 
character  line  Images.  Vertical  solid  lines  divide  cards  into 
input  fields.  Vertical  dashed  lines  in  input  fields  of  real 
numbers  specify  the  position  of  the  assumed  deci.mal  point  if 
the  decimal  point  is  omitted. 

4 . 1 Scene  Header 

Scene  header  [Fig.  ^.1]  is  a set  of  three  cards.  These 
cards  contain  scene  scale,  title  and  five  computational  tol- 
erances. These  three  cards  must  be  the  first  cards  of  every 


scene . 
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The  first  card  contains  delimiter  "1"  in  the  first  column. 
LIST  is  a diagnostics  flag.  A non-zero  value  will  cause 
listing  of  all  tables  after  the  scene  has  been  entered.  SCALE 
gives  the  number  of  straight-line  approximations  of  quadric 
curves  per  inch.  Note  that  small  value  will  give  the  drawn 
pictures  a coarse  appearance.  Large  value  will  require  large 
amount  of  memory.  Acceptable  value  of  SCALE  has  been  found 
to  be  15.0-20.0.  TITLE  is  a string  of  up  to  ^0  characters 
giving  the  name  of  the  scene. 

The  second  card  contains  five  computational  tolerances  : 
EPSl,  EPS2,  EPS3j  EPS4  and  EPS5.  Acceptable  values  of  these 
tolerances  have  been  found  to  be:  0.001,  0.0001,  0.00001, 
0.0001  and  0.001  respectively. 

The  third  card  is  a delimiter  card  with  "2"  in  the 
first  column. 

^ . 2 Surface  Equations 

Every  surface  equation  is  specified  by  a pair  of  cards 
[Fig.  4.2].  NS  is  the  number  of  a surface.  It  is  a positive 
integer  number  from  the  set  { 1 , 2 , . . . ,NSURFS } for  a scene  with 
with  NSURFS  surfaces.  The  surface  numbers  are  important 
because  each  surface  is  referred  to  by  its  NS  number  in  the 
remaining  sections  of  scene  specifications.  ITYPE  is  the  type 
of  the  surface:  1 is  component  surface,  2 is  auxiliary  sur- 
face. Auxiliary  surfaces  are  used,  for  example,  as  bounds  to 
eliminate  the  unwanted  sections  of  cones  or  hyperboloids  of 
two  sheets.  They  are  not  drawn  in  pictures.  The  equation  of 


Fig.  H.2  Surface  Equations 
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the  surface  is  specified  by: 

2 2 2 

a^x  +a2y  +a^xy+a^xz+agyz+a^x+agy+a^z+a^Q  = 0 

Mote  that  the  signs  of  the  coefficients  are  important.  Any 
point  outside  the  surface  must  yield  a positive  result  when 
substituted  to  the  above  equation.  Any  point  inside  the 
surface  must  yield  a negative  result. 

The  last  pair  of  surface  equation  cards  is  followed  by 
a delimiter  card  with  "3"  in  the  first  column. 

^ . 3 Surface  Bounds 

Surface  bounds  are  entered  in  format  shown  in  Fig.  ^.3- 
For  each  bounded  surface  there  is  at  least  one  card  giving  the 
number  MS  of  the  surface.  The  number  is  followed  by  up  to  15 
bound  values.  If  more  bound  values  are  give.n,  they  are 
entered  on  subsequent  cards  which  have  NS  omitted.  A list  of 
/ bounds  must  be  terminated  by  the  value  "9999".  If  a surface 

is  self-bounding,  such  as  a sphere,  its  bound  card  is  omitted. 
The  last  bound  card  is  followed  by  a delimiter  card  with  a "i*" 
in  column  one.  If  all  surfaces  in  a scene,  such  as  a set  of 
spheres,  are  self-bounded,  only  the  delimiter  card  is  entered. 

Each  bound  value  is  a non-zero  Integer.  The  absolute 
value  of  this  integer  is  the  surface  number  of  a surface  which 
bounds  surface  NS.  A positive  value  means  that  surface  NS  is 
j bounded  - that  is,  its  existence  is  limited  - to  the  space 

I 

I outside  the  bounding  surface.  A negative  value  means  that 

surface  NS  is  bounded  - its  existence  is  limited  - to  the  space 
j inside  the  bounding  surface.  The  "outside"  and  "inside"  space 


9 9 9 9 9 9 9 9 9 919  919  9 9 9 9|9  9 9 9 SN  9 9 9 
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of  a surface  is  determined  from  the  surface  equation  and 
therefore  the  signs  of  the  coefficients  of  the  equation  are 
Important.  A surface  may  have  several  sets  of  bounds  which 
are  separated  by  the  bound  value  ”9998".  Any  point  which 
satisfies  all  the  bounds  of  any  set  of  bounds  is  a point  on 
the  surface  of  NS.  This  indicates  that  the  bound  "9998" 
stands  for  logical  OR  operator  and  that  there  are  implied 
logical  AND  operators  between  bounds  in  each  set. 

4 . 4 Surface  Intersections 

There  is  one  card  as  shown  in  Fig.  4.4  for  each  inter- 
section of  two  surfaces.  The  last  intersection  card  is 
followed  by  a delimiter  card  with  a "5"  in  column  one.  If  a 
scene,  such  as  a set  of  non-intersecting  spheres,  does  not 
have  any  surface  Intersection,  only  the  delimiter  card  is 
entered. 

MSI  and  NS2  are  the  surface  numbers  of  the  two  inter- 
secting surfaces.  ISTYPE  gives  the  type  of  the  intersection 
as  follows; 

±1  two  planar  surfaces 

±2  one  planar  and  one  non-planar  surfaces 

±3  two  non-planar  surfaces,  however,  the  intersection 
is  planar 

±4  two  non-planar  surfaces 

>0  obtrusive  intersection,  that  is,  the  angle 

between  the  two  surfaces  Inside  the  object  is 
convex 


<0  intrusive  intersection,  that  is,  the  angle  between 
the  two  surfaces  inside  the  object  is  concave 
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If  either  of  the  two  intersecting  surfaces  is  non-planar, 
equation  of  a plane  which  intersects  the  intersection  must 
be  given.  The  equation  is  in  the  form: 
c^x  + + c^z  + c^^  = 0 

The  plane  is  referred  to  as  a cutting  plane  because  it  cuts 
the  intersection  and,  therefore,  determines  the  point  or 
points  on  the  intersection  where  the  program  will  begin  to 
trace  the  intersection. 

4 . 5 Real  Vertices 

There  is  one  card  for  each  real  vertex  or  a set  of  real 
vertices  as  shown  in  Fig.  4.5.  The  last  real  vertex  card  is 
followed  by  a delimiter  card  with  a "6"  in  column  one.  If  an 
object,  such  as  a cylinder,  does  not  have  any  real  vertices, 
only  the  delimiter  card  is  entered. 

M is  the  number  of  real  vertices  in  the  object  which  lie 
at  the  intersections  of  the  specificed  surfaces  MS^,...NS^^. 
There  must  be  at  least  three  surfaces  NS^,  NS2  and  NS^  to 
determine  a real  vertex  or  a set  of  real  vertices.  Maximum 
allowed  number  of  Intersecting  surfaces  is  15-  If  there  are 
less  than  15  Intersecting  surfaces,  the  last  one  must  be 
followed  by  "9999". 

Note  that  one  of  the  surfaces  NS^,  NS2  or  NS^  must  be 
planar.  If  none  is  available,  you  have  to  add  to  the  object 
an  auxiliary  planar  surface  which  Intersects  the  vertices  and 
specify  this  surface  as  one  of  the  first  three  surfaces  on 


the  vertex  card. 
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1 


i 

\ 
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^ . 6 Vantage  Points 

Each  vantage  point  Is  specified  by  a single  card  as  shown 
in  Fig.  4.6.  The  last  vantage  point  may  be  followed  by  an 
end-of-life  mark  or  by  a delimiter  card  with  a "9''  in  the 
first  column.  If  the  delimiter  card  is  entered  it  may  be 
followed  by  an  end-of-file  mark  or  by  the  cards  of  another  scene 
beginning  with  the  title  card. 

Each  vantage  point  card  contains  these  fields: 

LIST  is  a diagnostics  flag.  A non-zero  value  will  cause  the 
contents  of  all  tables  to  be  printed  after  the  vantage 
point  has  been  processed. 

IPTYPE  is  the  type  of  the  drawing  to  be  made  by  the  program: 

±1  all  lines  are  drawn 

±2  hidden  lines  are  dashed 

±3  only  visible  lines  are  drawn 

>0  orthographic  projection 

<0  perspective  projection 

PSCL  IS  the  picture  scale  of  the  drawing.  Ratio  of  the 

picture  scale  to  the  scene  scale  will  determine  the  size 
of  the  drawing. 

PHI,  THETA  and  PSI  are  the  azimuth,  elevation  and  twist  angles 
of  the  vantage  point. 

DE  is  the  distance  from  the  vantage  point  to  the  origin  of 
the  object  coordinate  system. 

DEE  is  the  distance  from  the  vantage  point  to  the  picture  plane. 


DE  and  DEE  are  used  for  perspective  projection  only.  They 
are  ignored  for  orthographic  projection  and,  therefore,  may 


be  omitted. 
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5.  DESIGN  OF  A QUADRAW  SCENE 

In  this  section  we  demonstrate  how  to  design  a QUADRAW 
scene  by  building  the  cup  shown  in  Fig.  5.1-3.  The  first 
part  of  the  data  file  consists  of  three  cards:  title  card 
with  scene  scale  of  15.0  and  a title  string;  tolerance  card 
with  five  tolerances  EPSl-5;  and  a delimeter  card. 

1 1 15.0  rUP  <NICHtEL  POTMESU) 

0.001  O.OCl  0.00001  0.0001  O.UOl 

2 _ 

Next,  we  have  to  specify  equations  of  surfaces,  surface 
numbers  [Fig.  5.1]  and  surface  types: 


Number 

Type 

Equation 

* 

Descript  ion 

1 

1 

2^2 

X + y - 16  = 

0 

Outer  cylinder 

2 

1 

2 2 

-x‘^  - y^  + 12. 

25  =0 

Inner  cylinder 

3 

1 

z - 6 = 0 

Top 

1 

z + 5.5  = 0 

Inner  bottom 

5 

]. 

-z  - 6 . 0 = 0 

Outer  bottom 

6 

1 

X 

t 

o 

II 

o 

Front  of  handle 

7 

1 

-X  - 0.5  = 0 

Back  of  handle 

8 

1 

2.56y^  + z^  - 

23.04y 

Outer  surface 

-z  + 36.09  = 0 

of  handle 

9 

1 

-iJ.Oy  - z*^  + 

36y 

Inner  surface 

+ z - 72.25  = 0 

of  handle 

10 

2 

y = 0 

Auxiliary  plane 
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Note  how  the  signs  of  coefficients  are  set  up.  Surfaces  1 
and  2,  for  example,  have  coefficients  with  opposite  signs 
since  the  space  outside  the  object  is  "outside"  of 
surface  1 and  "inside"  surface  2.  Surface  10  is  an  auxiliary 
transparent  surface  which  cuts  the  cup  in  two  halves  and  is 
used  to  eliminate  ambiguities  in  surf'ace  bounds.  We  enter 
these  cards  specifying  surface  equations: 


1 

1 

l.CO  1.00 

-16.00 

CUrSlDf  CYL 

2 

1 

o 

0 

1 

O 

0 

• 

1 

12.26 

INSIOE  CYL 

3 

1 

.1.00 

-6.0C 

TOP  Plr 

4 

1 

1.00 

5.60 

INSIOF  BOTTOM 

} 

1 

-1.00 

-6.  00 

OUTSIDE  BOTTOM 

6 

1 

1.00 

-0.50 

handle  80UN0 

7 

1 

-1.00 

-0.60 

HANDLE  BOUND 

a 

1 

z.st 

1.00 

-23.0'i 

-1.00 

36.09 

QUTSIDT  handle 

9 

1 

O 

0 

• 

1 

-1.00 

36.00 

1.00 

-72.26 

INSIOF  handle 

10 

2 

1.00 

0.00 

AUX.  PLANE 

Next,  the  surface  bounds  have  to  be  specified.  Surface 
bounds  are  used  to  specify  the  portions  of  a surface  which 
are  part  of  the  object.  For  example,  surface  number  3 - the 
top  of  the  cup  - is  a planar  surface  specified  by  z - 6 = 0. 
However,  only  the  portion  of  the  surface  which  lies  between 
the  outer  and  inner  cylinder  is  a part  of  the  object. 
Therefore,  surface  3 is  bounded  by  surfaces  1 and  2.  The 
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bounded  portion  of  surface  3 lies  "inside"  surfaces  1 and  2 
and  the  list  of  bounds,  therefore,  is  -1  -2  9999-  Bounds 
can  be  logically  grouped  by  the  special  bound  9998  which 
stands  for  logical  OR;  logical  AND  is  implied  between  two 
bounds;  9999  terminates  a list  of  bounds.  Surface  1 - 
the  outer  cylinder  - has  complicated  -bounds  because  it  is 
intersected  by  the  four  surfaces  of  the  handle.  It  is  also 
bounded  by  the  top  and  bottom  planes  and  an  auxiliary  surface 
number  10  is  used  to  limit  the  bounding  effects  of  surfaces  5 
and  6 to  the  right  half  of  the  cylinder  only.  We  obtain  a 
rather  complicated  list  of  bounds  for  surface  1: 

(-3  AMD  -5  AND  -10)  OR  (-3  AND  -5  AND  6 AND  10) 

OR  (-3  AND  -5  AND  7 AND  10)  OR  (-3  AND  -5  AND 
8 AND  10)  OR  (-3  and  -5  AND  9 AND  10) 

This  limits  surface  1 to  the  points  which  are  between 
the  top  and  bottom  planes  and  outside  the  two  handle  inter- 
sections. QUADRAW  limits  the  structure  of  bound  lists  to 
only  one  OR-AND  level  as  shown  above.  Similarly,  we  have  to 
carefully  check  all  the  other  surfaces,  except  the  auxiliary 
one,  and  for  this  cup  we  derive  the  following  bounds: 

1 -J  -5  -XO  9999  -3  -5  6 10  9998  -3  -i  7 10  <’998  -3 

-5  8 10  9999  -3  -5  9 10  9999 


2 

-3 

9 

9999 

3 

-1 

-2 

9999 

9 

2 

9999 

9 

-1 

9999 

6 

1 

10 

-8 

-9 

9999 

7 

1 

10 

-8 

-9 

9999 

e 

1 

10 

-6 

-7 

9999 

9 

1 

10 

-6 

-7 

9999 

9 
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Followlng  surface  bounds,  the  surface  intersections  are 
described.  For  example,  surfaces  1 and  3 intersect  in  an 
obtrusive  intersect  ion . The  intersection  curve  is  a circle; 
the  type  is  +2  and  a cutting  plane  is  x=0.  The  same  applies 
to  the  other  three  circular  intersections.  Notice  that  the 
intersection  of  surfaces  2 and  4 -has  type  -2  because  the 
intersection  is  intrusive  - the  angle  between  the  two  surfaces 
inside  the  object  is  greater  than  l80°.  Intersections  of 
the  handle  surfaces  and  the  outer  cylinder  do  not  need  any 


cutting  planes  since  they  are  terminated  by  real  vertices 


[Fig. 

5. 

2j. 

The  specified  surface  intersections  are: 

1 

i 

2 

l.OO 

1 

5 

2 

1.00 

z 

3 

•j 

1.00 

2 

-2 

1.00 

1 

6 

-2 

1 

7 

-2 

1 

P 

-2 

1 

9 

-2 

6 

9 

2 

6 

9 

2 

7 

P 

2 

7 

5 

9 

2 

The 

last  part  of  scene  design  is  to  specify  the  real 

vertices 

of 

each  object.  A real  vertex  is  a point  on  the  sur 

face 

of 

the 

object  where  three  or  more  surfaces  Intersect. 

There 

are  8 

real  vertices  in  the  cup.  They  form  4 sets  of 

real 

vertices  - each  set  has  2 real  vertices.  Each  set  lies 

at  the  intersection  of  the  outer  cylinder  - surface  number  1 - 

and  two 

handle  surfaces  [Fig.  5*2].  We  enter  these  cards: 

2 

a 

6 

1 9999 

Z 

a 

7 

1 9999 

2 

9 

6 

1 9999 

2 

9 

7 

1 9999 

6 
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Following  the  scene  data  we  add  several  vantage  point 
specifications  and  terminate  the  file  by  a "9"  delimiter 
card : 


3 

20. CO 

•i  5.00 

‘>5,00 

22.  50 

3 

20.00 

22.00 

22.00 

0.00 

3 

20.00 

70. LO 

10.00 

30.00 

3 

20.00 

20.00 

60.00 

-100.00 

3 

20.00 

0.00 

22.50 

22.  50 

2 

20.00 

C.OO 

22.50 

30.00 

3 

20.00 

100.00 

-30.00 

80.30 

2 

20.00 

‘.5.00 

75.50 

30.00 

9 
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Appendlx  A 

CDC  6600  Version  of  QUADRAW 


QUADRAW  has  been  run  on  a CDC  6600  computer  under  the 
KRONOS  2.1  or  NOS  1.1  operating  system  and  compiled  by  an 
FTN  4.3  compiler. 


1 


A . 1 FORTRAN  Source  Text  Library 

The  FORTRAN  source  text  of  the  program  Is  stored  as  a 
library  file.  It  is  maintained  by  the  CDC  text  editor  UPDATE. 
The  text  library  contains  CDC,  Univac  and  IBM  versions  of  the 
program.  Versions  for  other  computers  can  be  added  to  this 
library.  UPDATE  can  write  on  file  COMPILE  either  the  text  of 
the  CDC,  Univac  or  IBM  version.  The  CDC  version  is  compiled 
by  FTN  from  COMPILE  into  a binary  file.  The  Univac  or  IBM 
version  is  copied  from  COMPILE  to  a magnetic  tape  and  trans- 
ported to  a Univac  or  IBM  computer. 

The  computer-dependent  parts  of  the  master  text  file  are 
written  as  follows: 

FORTRAN  source  text  common  to  all  versions 
»IF  DEF,CDC 

FORTRAN  source  text  for  CDC  version  only 
»ENDIF 

*IF  DEF, UNIVAC 

FORTRAN  source  text  for  Univac  version  only 

»ENDIF 
»IF  DEF, IBM 

FORTRAN  source  text  for  IBM  version  only 
»ENDIF 


FORTRAN  source  text  common  to  all  versions 
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These  conditional  IP-ENDIF  blocks  are  Inserted  whenever  there 
are  computer-dependent  sections  of  QUADRAV/.  Text  editor 
UPDATE  will  write  the  conditional  blocks  on  file  COMPILE  as  a 
part  of  the  text  to  be  compiled  whenever  the  conditional  value 
of  the  block,  i.e.,  CDC  , UNIVAC  or  IBM,  is  defined. 

The  QUADRAW  source  text  is  also  maintained  on  an  Adage 
computer  as  a back-up  copy  [5].  It  can  be  written  on  magnetic 
tape  in  various  formats  [6]. 

A. 2 QUADRAW  Execution  on  CDC  6600 

A QUADRAW  magnetic  tape  contains  two  UPDATE  libraries 
written  as  two  tape  files.  The  first  file  is  the  master 
source  text  of  the  program  QUADRAW.  The  second  file  is  the 
data  library  containing  specifications  of  several  QUADRAW 
scenes.  The  magnetic  tape  is  assigned  to  a Job  by  the 
following  statement : 

LABEL(TAPEIN,VSH=Txx,MT,D=HY,F=SI,LB=KU,PO=R) 
where  the  parameters  are: 


TAPEIN 

Local  file  name 

VSN=Txx 

Tape  VSW  number 

MT 

7-track 

D=HY 

800-bpi 

F=SI 

Format : Scope  Internal 

LB=KU 

Labels:  Kronos  Unlabelled 

II 

o 

0. 

Read  only 

a)  To  copy  the  first  file,  update  It,  compile  it  and 
catalog  the  new  UPDATE  library  and  binary  file: 


LABEL (TAPEIN,  . ) 

REWIND (TAPEIN) 

COPYBF(TAPEIN ,OLDPL, 1 ) 

PURGE (QLIB/NA) 

DEFINE ( NEWPL=QLIB ) 

UPDATE ( F , P=OLDPL , N=MEWPL ) 

FIN(I=C0MPILE,L=0) 

PURGE  (BQDRAV;/NA) 

DEFINE (BQDRAW) 

COPYEI ( LGO , BQDRAW , X ) 
end-of-record 
*DEFINE  CDC 
*IDENT  PLOTS 
»DELETE  PL0TS3L.1 

CALL  PLOTS (LIMIT plotter  output  ident ' ) 
end-of- information 


This  run  will  produce  two  permanent  files: 

QLIB  master  text  library  of  QUADRAW 

BQDRAW  binary  file  of  WUADRAW 

b)  To  copy  and  update  the  second  file  - QUADRAW  scenes: 


LABEL  (TAPEIN, ) 

REWIND (TAPEIN) 

SKIPP(TAPEIN,1) 

COPYBF (TAPEIN , OLDPL , 1 ) 

PURGE (QSLIB/NA) 

DEFINE ( NEWPL=QSLIB ) 

UPDATE ( P=OLDPL , N=NEWPL ) 
end-of-record 
* IDENT  NEW 
end-of- in format  ion 

This  run  will  produce  a permanent  file  QSLIB  contain 
Ing  several  DECKS.  Each  DECK  contains  the  specif lea 
tions  of  a QUADRAW  scene  and  several  vantage  points. 
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c)  To  execute  QUADRAW ; input  - FORTRAN  unit  5 - is  read 

from  file  COMPILE:  output  - FORTRAN  unit  6 - is  written 
on  file  OUTPUT: 

NAME,T100,CM15000-  QUADRAW 

ACCOUNT  ( . ... ) 

CHARGE  ( . . . , ) 

REV;iND(  OUTPUT) 

ATTACH ( OLDPL=QSLIB ) 

UPDATE (P=OLDPL,D) 

ATTACH (BQDRAW) 

LOAD(BQDRAW) 

EXECUTE ( QDRAW , COMPILE ) 
end-of-record 
*COHPILE  CUP 
*COMPILE  FIG2 
end-of-in format  ion 

The  COMPILE  file  must  have  80  characters  of  data  per 
line;  therefore,  use  the  D option  on  the  UPDATE  card 
as  shown  above.  Two  or  more  decks  of  QUADRAW  scenes 
from  file  QSLIB  may  be  combined  by  UPDATE  on  file 
COMPILE  and  run  together. 


d)  The  Univac  or  IBM  version  of  QUADRAW  are  generated 
from  the  master  source  text  on  a CDC  6600  by  the 
following  run; 


LABEL  (TAPEOUT, ) 

REWIND (TAPEOUT) 

ATTACH (OLDPL=QLIB) 
UPDATE(F,P=OLDPL) 

C0PYCF( COMPILE , TAPEOUT , 1 ) 
ATTACH (0LDPL=QSLI3) 

REWIND (COMPILE) 
UPDATE(P=OLDPL,D) 

C0PYCR( COMPILE .TAPEOUT , 10 ) 
REWIND (TAPEOUT) 
end-of-record 

»IDENT  UNIVAC  or  »IDENT  IBM 

end-of-record 

♦COMPILE  APOLLO 
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*WEOR 

^COMPILE  CUP 
*WEOR 

*COMPII.E  SPHR8 
*WEOR 


end-of-inf ormat ion 

Change  the  density,  parity  and  7/9  track  parameters 
in  the  above  LABEL  statement  as  required  by  the 
Uni vac  or  IBM  computer. 

A. 3 SCOOP  Plotting  Package 

QUADRAW  uses  the  following  subroutines  of  the  SCOOP 
plotting  package: 

PLOTS(LIMIT, 'your  identification  string') 
PLOT(X,Y,IPEN) 

PLOTS  opens  plotting  output  up  to  LIMIT  inches  long  and 
labels  it  with  'your  identification  string'. 

PLOT  moves  the  pen  from  the  current  location  to  location 
(X,Y).  IPEN  can  have  the  following  values: 


2 

solid 

line (pen 

down) 

and 

do  not  reset  plot 

origin 

3 

skip 

(pen 

up 

) 

and 

do  not  reset  plot 

origin 

-2 

solid 

line (pen 

down ) 

and 

reset  plot  origin 

to  (X,Y) 

-3 

skip 

(pen 

up 

) 

and 

reset  plot  origin 

to  (X,Y) 

999 

skip 

(pen 

up 

) 

and 

close  plotting  output 

A. ^ Timing  Routine 

MSTIME  calls  PTM  library  subroutine  SECOND(TI.ME ) which 


returns  the  time  since  the  beginning  of  program  execution  in 
seconds  in  floating  point  variable  TIME. 


r — ^ 
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Notes 

The  current  CDC  implementation  of  QUADRAW  has  these 
features: 

' a)  Each  picture  is  drawn  on  one  12  x 12  inch  sheet 

of  paper. 

b)  PLOTS  opens  output  on  green  graph  paper. 

PLOTSBL  with  the  same  arguments  as  PLOTS  opens 
output  on  white  blank  paper. 

c)  LIMIT,  i.e.,  the  maximum  total  length  of  the  paper, 
is  set  to  600  inches. 

51 
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Appendlx  B 

Unlvac  1110  Version  of  QUADRAW 

QUADRAW  has  been  run  on  a Univac  1110  computer,  in 
demand  mode  under  EXEC -8. 

There  are  two  types  of  magnetic  tape  files  from  which 
QUADRAV/  can  be  run  on  a Univac  computer: 

a)  Univac  program  file  written  in  COPY  format  con- 
taining symbolic  and  relocatable  elements  of  all 
QUADRAW  routines  and  symbolic  elements  of  all 
QUADRAW  scenes. 

b)  File  written  in  BCD  code  on  a CDC  6600  computer. 

This  file  is  read  and  converted  by  a Univac  as 
described  in  [6], 

3. 1 A QUADRAW  Run 

The  following  EXEC-8  statements  illustrate  how  to  run 
QUADRAW  on  a Univac  1110  computer  from  a demand  terminal. 
Assume  that  we  have  a 9-track  QUADRAW  tape  number  Txxxxx 
v;ritten  in  Univac  COPY  format.  We  assign  an  unregistered 
mass  storage  file  called  QDRAW  and  copy  the  tape  file. 

@ASG,UP  QDRAW. 

@MSG,W  MAY  I USE  TAPE  ’Txxxxx'  9-TRK  FOR  5 MIN 
@ASG,T  TAPEIN . ,T jTxxxxx 

SCOPYjG  TAPEIN. , QDRAW. 

5FREE  TAPEIN. 

Pile  QDRAW  contains  all  symbolic  and  relocatable  elements 
of  QDRAW.  We  may  list  the  directory  of  QDRAW: 


gPRT,T 


QDRAW. 
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If  It  is  necessary  to  edit  or  compile  any  element  use 
the  editor  and  compiler: 

SED  QDRAW. eltname ,QDRAW. eltname 

@FOR  QDRAW. eltname , QDRAW. eltname 

When  all  relocatable  elements  are  updated  we  create  an 
absolute  element: 

3 PREP  QDRAW. 

@MAP  QDRAW . QDRAW , QDRAW . QA3S 
IN  QDRAW. 

LIB  GRAIL$*LI3. 

END 

An  absolute  element  called  QABS  was  added  to  file  QDRAW. 

This  element  can  be  executed  after  the  plotting  package  was 
directed  to  plot  its  output  on  a digital  plotter: 

3GRAIL$*LIB . PLOTMODE ,C 

§XQT  QDRAW. QABS 

QUADRAW  is  now  waiting  for  the  first  scene  to  be  entered. 

V/e  may  enter  the  CUP  and  APOLLO  scenes  by: 

3ADD  QDRAW. CUP 

3ADD  QDRAW. APOLLO 

SEOF 

The  execution  of  the  program  was  terminated.  The  file 
QDRAW  may  be  saved  back  on  the  same  tape  as  the  second  tape 
file  by  com.mands: 

1 

t 

I 

I 

i 
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SASG,T  TAPEOUT. ,T,TxxxxxR 
@TSKIP  TAPEOUT. 

3COPY,G?'I  QDRAW,  , TAPEOUT. 

3FREE  TAPEOUT. 

B. 2 GRAIL  Plotting  Package 

QUADRAW  uses  the  following  subroutines  of  the  GRAIL 

plotting  system  at  the  Univac  1110  installation: 

SETPLT  (300.0,11.0,0) 

ORIGIN (XORG , YORG ,M0DE ) 

SCRIBE(X,Y) 

SLEW(X,Y) 

ENDPLT 

SETPLT  opens  plotting  output  and  establishes  300  x 11 
inch  maximum  plotting  area.  ORIGIN  positions  the  pen  to 
the  center  of  a new  11  x 11  inch  picture  area.  SCRIBE  draws 
a line  from  the  current  position  of  the  pen  to  position  X,Y. 

SLEW  moves  the  pen  (in  the  up  position)  to  position  X,Y. 

ENPLT  closes  the  plotting  output. 

B. 3 Timing  Routine 

MSTIME  calls  the  system  subroutine  CPTI;'E(M)  -which  returns 
the  elapsed  execution  time  of  the  run  in  integer  milliseconds 
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Appendix  C 

IBM  System/360-370  Version  of  QUADRAW 

QUADRAW  has  been  run  on  an  IBM  36O/67  computer  using 
the  Alpha  time-sharing  system. 

The  IBM  version  of  QUADRAV/  is  written  on  magnetic 
tape  by  a CDC  computer  in  BCD  code.  It  is  read  and  listed 
by  an  IBM  computer  with  the  following  JCL  statements: 


//QDRAW 

JOB.  . 

//STEPl 

EXEC 

PGM  = IE3PTPCH,REaiOM=32K 

.'■/SYSPRIXT 

DD 

SYS0UT=A 

//SYSIM 

DD 

DUMMY 

//SYSUTl 

DD 

UNIT=( TAPE, 1, DEFER) 

//SYSUT2 

DD 

V0L= ( , RETAIN , 1 ,SER=xxxxxx ) 

DSN= QDRAW, 

LABEL=(1,NL) 

DISP= (OLD, KEEP, KEEP) 

DC3= ( RECFM=FB , LRECL=80 , 3LKSIZE=80 
3UFN0=2 ,DEN=2 ,TRTCH=ET ) 
sysouT=A 

//SYSIM 

DD 

* 

/* 

To  punch  the 

' file 

on  cards  use  this  SYSUT2  statement 

//SYSUT2 

DD 

sysouT=3 

CALCOMP  Plot 

ting 

Package 

QUADRAV/  uses  the  folio  wing  subrout  ines  of  the  CALCOMP 


plotting  package: 

PLOTS (0,0,0) 
PLOT(X,Y,IPEN) 


PLOTS  opens  plotting  output.  The  plotted  paper  is 
labelled  at  both  ends  with  the  Job  name. 

PLOT  moves  the  pen  from  the  current  location  to  location 
(X,Y).  IPEN  can  have  the  following  values: 
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2 solid  line  (pen  down)  and  do  not  reset  plot  origin 

3 skip  (pen  up  ) and  do  not  reset  plot  origin 

-2  solid  line  (pen  down)  and  reset  plot  origin  to  (X,Y) 

-3  skip  (pen  up  ) and  reset  plot  origin  to  (X,Y) 

999  skip  (pen  up  ) and  close  plotting  output 

C . 2 Timing  Routine 

MSTIME  calls  system  routine  TIMEH(M)  which  returns  the 
time  since  the  beginning  of  program  execution  in  tenths  of 
seconds  in  integer  variable  M. 

Notes 

The  current  IBM  implementation  has  these  features: 

a)  Each  picture  is  drawn  within  an  3 x 3 inch 
area  of  paper. 

b)  The  plotter  uses  23-inch-wide  paper.  Therefore, 
three  pictures  are  drawn  along  the  width  of  the 
paper . 
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Appendix  D 
QUADRAW  SCENES 


A set  of  ten  scenes  has  been  used  to  test  the  program. 

Each  scene  description  is  followed  by  several  vantage  point 
specifications  which  are  terminated  by  a "9"  delimiter  card. 
Several  scenes  can  be  combined  for  a single  program  run. 

The  scenes  for  which  descriptions  exist  are: 

APCLLO  - a model  of  the  Apollo  spacecraft  [Fig.  D.1-2]. 

CUP  _ a model  of  a cup  [Fig.  5.1-3]- 

SPHRS  - a scene  made  of  8 spheres  which  are  cut  by 

two  planes  [Fig.  1.2  and  D.3]. 

CYLIN4  - a scene  made  of  4 cylinders  [Fig.  D.4]. 

FIGl  - a sphere  intersecting  with  a cylinder 

[Fig.  D.  5]. 

FIG2  - two  cone  parts  intersecting  with  a cylinder 

[Fig.  D . 6 J . 


ELPHYPl  - four  elliptic  hyperboloids  of  one  sheet 
[Fig.  D .7 ] . 

ELPHYP2  - four  elliptic  hyperboloids  of  two  sheets  [Fig.D.8].  j 
CYLIN4P  - four  cylinders  with  perspective-projection 
vantage  points  [Fig.  D.9]. 

SPHERES  - a scene  made  of  eight  non-intersecting 

spheres  [Fig.  D.IO].  : 

I : 

I' 
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1. 

CLOSER 

2.3.11 

2. 

CUT 

2.3.1 

3. 

DELEV 

2.3.6 

il. 

DERIV 

2.5.9 

5. 

DET3 

2.5.10 

6. 

DIRECT 

2.3.5 

7. 

DRAW 

2.2.10 

8. 

ELX21 

2.5.7 

9. 

ELY2 

2.5.8 

10. 

EQlll 

2.5.1 

11. 

EQ211 

2.5.1 

12. 

EQ22 

2.5.2 

13. 

EQ221 

2.5.1 

14. 

EXTRM 

2.3.8 

15. 

FINAL 

2.2.9 

16. 

FSTPT 

2.3.7 

17. 

FULL 

2.4.6 

18. 

GTOll 

2.4.1 

19. 

GT012 

2.4.1 

20. 

GT013 

2.4.1 

21. 

GT031 

2.4.1 

22. 

GT041 

2.4.1 

23. 

GT051 

2.4.1 

24. 

GT061 

2.4.1 

25. 

GT069 

2.4.1 

26. 

GT071 

2.4.1 

27. 

GT072 

2.4.1 

28. 

GT082 

2.4.1 

29. 

GT091 

2.4.1 

30. 

GT151 

2.4.2 

31. 

IBND 

2.3.2 

32. 

IBND2 

2.3.3 

33. 

INITO 

2.3.12 
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3^. 

INPOB 

2.2.1 

35. 

IMPVW 

2.2.4 

36. 

INSIDE 

2.2.8 

37. 

INTBL9 

2.4.8 

38. 

INVIS 

2.3.9 

39. 

IPNT 

2.5.12 

40. 

MSTIME 

2.7.1 

41. 

NEWT 

2.5.5 

42. 

POLAR 

2.5.11 

^3. 

PWRT 

2.5.6 

44. 

POLRT 

2.5.6 

45. 

QDRAW 

2.1.1  (main  program) 

46. 

REDGE 

2.2.6 

47. 

ROOT2 

2.5.3 

48. 

ROTPT 

2.3.4 

49. 

RVERT 

2.2.2 

50. 

SINT 

2.2.3 

51. 

SORTl 

2.3.13 

52. 

STOll 

2.4.3 

53. 

ST012 

2.4.3 

54. 

ST021 

2.4.3 

55. 

ST031 

2.4.3 

56. 

ST032 

2.4.3 

57. 

ST041 

2.4.3 

58. 

ST042 

2.4.3 

59. 

ST051 

2.4.3 

60. 

ST061 

2.4.3 

61. 

ST071 

2.4.3 

62. 

ST072 

2.4.3 

63. 

ST073 

2.4.3 

64. 

ST081 

2.4.3 

65. 

ST085 

2.4.3 

66. 

ST151 

2.4.4 

67. 

ST152 

2.4.5 

68. 

SUBSl 

2.5.4 
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69. 

TABLE 
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70. 

TRMSF 

2.2.5 

71. 

VEDGE 

2.2.7 

72. 

WDUM? 
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